Transient quantum evolution of 2D electrons under photoexcitation of a deep center 
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We have considered the ballistic propagation of the 2D electron Wigner distribution, which is 
excited by an ultrashort optical pulse from a short-range impurity into the first quantized subband of 
a selectively-doped heterostructure with high mobility. Transient ionization of a deep local state into 
a continuum conduction c-band state is described. Since the quantum nature of the photoexcitation, 
the Wigner distribution over 2D plane appears to be an alternating-sign function. Due to a negative 
contribution to the Wigner function, the mean values (concentration, energy, and flow) demonstrate 
an oscillating transient evolution in contrast to the diffusive classical regime of propagation. 

PACS numbers: 05.30.-d; 73.20.-r; 78.47,+p 



I. INTRODUCTION 

In recent decades, intensive efforts were paid in or- 
der to study the quantum coherent properties of differ- 
ent physical systems [lj. During the development of the 
ultrafast spectroscopy of bulk semiconductors and het- 
erostructures Q , both coherent oscillations between cou- 
pled states and different relaxation processes have been 
investigated (see references in Q and [I])- Some quan- 
tum peculiarities, e.g. in the transport of mesoscopic de- 
vices [j| or in the dynamics of electron excited at metallic 
surfaces [f|, were discussed recently but, to the best of 
our knowledge, the coherent dynamics of a free quasi- 
particle, which propagates over continuum states, was 
not measured directly in any solid state system. The 
quantum response, such as the formation of quasiparti- 
cles in different systems 0, H, H, 03 > has been observed 
for sub-picosecond stage of evolution. Under theoretical 
consideration of such kind of measurements (e.g., see 11 1 
and references therein), one can model the photogener- 
ation process using a simple initial condition describing 
the creation of carriers during a femtosecond temporal in- 
terval. At the same time, in the case of photoexcitation 
of carriers with low concentration and with energy values 
below the optical phonon energy, the dynamical regime of 
the response appears to be valid up to nanosecond time 
interval. It is because both the fast relaxation, due to op- 
tical phonon emission, and the carrier-carrier interaction 
are suppressed. Thus, a possibility is to study the quan- 
tum nature of the ballistic transient evolution, caused by 
the non-classical character of photoexcitation 0, E3] • 

Modern high-mobility heterostructures are character- 
ized by a momentum relaxation time correspondent to 
the subnanosecond scale at low temperatures So 
that the mean free path appears to be macroscopic (>100 
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/ini if electron energy is about few meV) . The photoex- 
citation of a single deep impurity under a laser pumping 
focused up to submicron scale [14| can be carried out in 
a non-doped heterostructure with a low surface concen- 
tration of centers (deep centers in bulk GaAs are under 
consideration since the starting of 70s [IB]). Below we 
consider the transient photoexcitation of electrons from 
a deep impurity level and the quantum ballistic evolution 
of the Wigner distribution in the 2D plane over submil- 
limeter distances during nanosecond time interval. 

In contrast to transitions between local states, when 
the Rabi effect (oscillations of population versus pump- 
ing intensity [1, [ll]) takes place, the level population of 
a local state under ultrafast photoionization decreases 
monotonically with the pumping intensity because the 
excited electron appears to be delocalized over the con- 
duction c-band. Another peculiarity of the process under 
consideration is the quantum character of the transient 
evolution. Due to this, the concentration distribution, 
which decreases from the center, involves an oscillating 
contribution and regions of a negative Wigner distribu- 
tion take place. Such a distribution should be considered 
with the use of the quantum kinetic equation, written in 
the Wigner representation, due to the following reasons: 
(a) the energy conservation law is not valid during the 
photoexcitation process and (b) there is no momentum 
restrictions on the excited distribution due to the short- 
range impurity state involved in the phototransition. 

In this paper, we restrict ourself to the local time ap- 
proximation which corresponds to the photoionization 
above the c-band edge (AuT ex > 1, where Aw is the 
detuning frequency and r ex is the duration of photoexci- 
tation), when only the point (b) is essential. Due to this 
reason, the mean values (concentration, energy, and flow) 
show an oscillating behavior in contrast to the diffusive 
classical regime. Moreover, although the concentration 
and energy distributions are positive-definite functions, 
the flow distribution appears to be an alternating-sign 
one, i.e. a flow may be directed opposite to a concen- 
tration gradient. The peculiarities discussed can be veri- 
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fled by the use of optical methods or scanning tunneling 
microscopy, if the measurements can be performed with 
submicron and subnanosecond resolutions. 

The present work is organized as follows. The photoex- 
citation process, including the evolution of the deep cen- 
ter population and the transient Wigner distribution over 
c-band, is described in Sec. II. Section III presents tem- 
poral dependencies of the above-introduced functions. 
The transient dynamics of the mean values is described 
in Sec. IV. A list of the assumptions used and the discus- 
sion of the methods for experimental verification of the 
peculiarities discussed are given in the concluding sec- 
tion. Appendix contains the description of the classical 
regime of transient evolution. 



e.p = p 2 /m with the effective mass of c-band, m, the 
form- factor w t introduced by the relation E t — Ew f , and 
the detuning energy TlAlo, which takes into account the 
subband quantization. 

We are using the initial conditions f PlP2 t^-oo = 
and nt—>-oo = 1, which correspond to the single-electron 
population of the spin-degenerated level, so that the nor- 
malization condition takes the form: n t + 2 / p , p t = 1. 
Evolution of the local state population is governed by the 
integro-differential equation 

dn t ( eEv cv \ 2 f* 
x^|/ p | 2 cos(|-A w )(t-O = 0, (4) 



II. ULTRAFAST PHOTOEXCITATION 

Under photoexcitation of electrons, transitions from a 
deep local level into the first subband of c-band is de- 
scribed by the density matrix of the j-state, ~Pj t . Per- 
forming the averaging over the period of the radiation 
Et exp(—iujt)+c.c. one obtains the quantum kinetic equa- 
tion 113: 



which is obtained from Eqs. (1, 2). Instead of Eq. 
(3), one can describe the transient evolution of the o 
band distribution through the Wigner function, / Pi qt = 
/p-ftq/2, P +aq/2t which is governed by the equation 
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Here hj is the Hamiltonian of the 2D state in the c-band 
(j = c) or of the state at short-range centre (j = h), 
Sjf—t = exp[—ihj(t' — t)/h] is the evolution operator of 
the j-th state, and (y)jj' is the velocity matrix element 
for j «-> j' transitions. For the case of a deep center con- 
nected to the valence w-band we use in (2) the inter- 
band matrix element of velocity, v cv , multiplied by the 
overlap integral between the plane wave of momentum 
p and the local state, I p = (p\h). The evolution of the 
distribution function over c-band, f PlP2 t = (PifPctlPa)) 
is governed by the equation: 



d/pip 



2 l 



... + h^ Spi £ P 2 )-^PiP 2 * — ^PiPa*' 



G 



PlP2* 



eEv c . 



I Pl I* 2 w t / dt'wf 

J — GO 



Xe i(e p /h-Au,)(t'-t) rhi + ( 



c.c, pi 



Pa) 



(3) 



with the right-hand side dependent on the population 
of the local state, n t = (h\p ht \h). The generation rate, 
G PlP2 t, is determined through the 2D kinetic energy 



with the velocity v = p/m. Similarly, the generation 
rate, G p , qt = G p _ h q/ 2 , P +h< 1 /2t, is transformed into 
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and the right-hand side of Eq. (5) is determined through 
the evolution of fit- The solution of (5) takes the form 
/ P ,qt = JL^ dt' exp[-zq • v(t - t')]G Pitlt >, so the descrip- 
tion of the transient evolution is reduced to the calcu- 
lation of a multiple integral and to the solution of Eq. 
(4). 

A simplified consideration of the problem is possible 
under the condition AuiT ex > 1 (photoionization into a 
high-energy state of c-band) when Wynv in Eqs. (4) and 
(6) can be replaced by w t n t due to the fast oscillating 
factors (the local time approximation). The integration 
over dt' in Eq. (6), which is performed with an infinites- 
imal damping factor in the exponent, <5 ^ +0, gives 



(-rp.qt ~ \ % ; j l p +hq/2 1 P -hc l /2W t n t 
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Here the kinetic energy of the inhomogeneous system 
(q 7^ 0) is not conserved during the photogeneration pro- 
cess, even for the long r ex case, due to the violation of 
the momentum conservation law. Using Eq. (7) and per- 
forming the Fourier transformation of the distribution 
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FIG. 1: Transient evolution of the population of short-range 
level under the dimensionless pumpings yr ex /2 =0.33, 1, 3 
and 9. 



function / p . q t one obtains: 
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,(8) 



where we have introduced the time-dependent coordi- 
nate, x f _e = x v(i — i'). The function / PiX t satis- 
fies the conditions /_ p ,- x4 = / p , X i and /_ p , xt = / p _ xt , 
which are verified by the exchange q — * — q in Eq. (8). 

Within the local time approximation, Eq. (4) takes 
the form: 



7 = Air 
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V flUJ 
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where 7 stands for the photoionization decrement. The 
analytical solution of Eq. (9), 



n t = 1 - 7 



y cft'i4 ex P ^-7 y 



(10) 



describes the transient population of the level [17[. Thus, 
we have obtained the Wigner distribution (8), and the 
population (10) written in the integral forms which con- 
tain the overlap integral I p and the form-factor Wt- 



III. SHORT-RANGE CASE 

To calculate the distribution (8) and the population 
(10) we use below the overlap integral for the short-range 
local state with the characteristic size l a , so that I p ~ 
l a /L if p < h/l a and I p ~ if p > h/l ; here L is the 
normalization length. Within the above assumption, the 
decrement of photoionization in Eq. (9) takes the form: 
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where p S2 d ^ s l ne density of states, is the heavy 
hole mass, and the level coupling energy, AE, is ex- 
pressed through l a according to AE ~ {h/l ) 2 /2m/,. The 
temporal dependencies of n t under different pumping 
level, which is determined by the dimensionless parame- 
ter 7 r ex /2, are shown in Fig. 1 for the Gaussian form- 
factor wt = exp[— 2(t/T ex ) 2 }. The complete ionization of 
the center appears under the condition "~)T ex /2 ~2 and, 
when 7 increases, the photoionization takes place dur- 
ing the front of pulse. The full ionization regime takes 
place under a pulse energy ~ 0.2/xJ focused on an area 
~ 100/xm; this estimate is performed for the GaAs pa- 
rameters and does not depend on the pulse duration. 

Next, we turn to the description of the photoexcited 
electron distribution given by Eq. (8) and dependent 
on time, |p|, |x|, and the angle pTx. We consider the 
long-duration excitation case (the dynamic regime of re- 
sponse takes place for the nanosecond time scale) and 
demonstrate that the Wigner distribution / PjX t is not a 
positive-definite function. We calculate below the dis- 
tribution at the maximal pumping, t — 0, for the cases 

p||x and p_Lx with the use of the notations fp} x and fpj^, 
respectively. Performing in Eq. (8) the integration over 
q_Lp one obtains the distributions fp}h ± as follows: 
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where 8[z] is the Heaviside step function, v — |v|, and 
PAu, = \l2rnhAui is the characteristic momentum. The 
integrals over t and p\ can be factorized for the slow 
electron case, p ~ 0; moreover, a non-zero contribution 
appears from the first addendum only. The distributions 
for the II and _L orientations are coincident, 



Jp=0,x 



2 7 



dtw t n t Jq 



2pAuiX 



(13) 



and the coordinate dependence is given by the zero- 
order Bessel function, Jq(z), which has an alternating- 
sign value and decreases as a square root. 
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The distribution functions (12) depend on the dimen- 
sionlcss momentum and coordinate, p/pauj and 
the detuning parameter, AwT ra , and the pumping in- 
tensity, 7. Using (10) with the dimcnsionless pumping 
jr ex /2 = 3 and performing the numerical integration in 
Eq. (12) we plot the Wigner distribution for Awr ex = 10 
as it is shown in Fig. 2. One can see a non-monotonically 
dependency on p/p Au> and xp^jh with pronounced neg- 
ative contributions. Both longitudinal and transverse 
cases show a fast decrease with dimensionless momen- 
tum, whereas oscillations slowly decrease with dimen- 
sionless coordinate due to the spread of the distribution 
under propagation. 



IV. MEAN VALUES 

The transient dynamics of the Wigner distribution un- 
der consideration can be verified by the treatment of 
spatio-temporal dependencies of the mean values (con- 
centration, energy, and flow, n xt , £ x t, and i xt ) given by 
the standard formulas: 



n x t 
ixt 



dp 



{2ttT%Y 



/pxi- 



(14) 



Below we analyze the spatio-temporal evolution of (14) 
using the distribution (8). Performing the integration 
over the variable p ± one obtains the concentration 
and energy distributions, which are isotropic over the x- 
plane: 
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where V means the principal value of the integral, Ji(z) 
is the first-order Bcsscl function, vau = PAu/m, and y = 
s/HAlo is the dimensionless energy. The distributions 
(15) depend on the pumping intensity through n t , given 
by Eq. (10), and on the dimensionless coordinate and 
time, xpAuj/h and t/r ex . 
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FIG. 2: Longitudinal (a) and transverse (b) Wigner distri- 
butions at maximal pumping (t = 0) versus dimensionless 
momentum and coordinate, p/pAui and xpAu/fi- 



Before numerical calculations, we consider the asymp- 
totes of n xt (similar formulas can be written for £ xt 
and the flow distribution) for the case t » r ex and 
x >• h/pAw Using the asymptotic expansion of the 
Bessel function for large arguments and performing the 
integrations over t' and q, one obtains the explicit ex- 
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where TV = dtw 2 n t /T ex , Iau = K/paw, and we 
have introduced the dimensionless forms z± = (x ± 
VAwt)pAu}/2hv / Aujt and z y± = z± ± (y/y - l)V Atut. 
Thus, the period of o scilla tions of the concentration dis- 
tribution, (?i/2pAw)v / Awt, does not depend on Aw and 
increases as \ft. 

The spatio-temporal dependency of the concentration 
is shown in Figs. 3. Fig. 3(a) shows the 3D graph 
near the peak at xpAw/h — for Au)T ex = 5 and 
jT ex /2 — 1. As can be seen, concentration falls quickly 
for small dimensionless coordinate and oscillates for big- 
ger xpAu/h values. At the same time maximum posi- 
tion in dimensionless time is shifted following the clas- 
sical velocity as shown in Fig. 3(6) where the contour 
plot together with the line corresponding to the classi- 
cal velocity is presented. Concentration is normalized by 
no = "fT ex m/\Lo/Tth, which is equal to 2.82 x 10 11 cm~ 2 
for %Au = 5 meV, and -fT ex /2 = 1. This value corre- 
sponds to an excited electron localized over an area of 
the order of (H/palu) 2 - 

Fig. 4 shows the behavior of the energy distribution 
vs. dimensionless position and time. Energy distribution 
has been normalized by E — n n hAui, which corresponds 
to 2.25 x 10~ 3 erg/cm 2 for the same values of HAlu and 
"fr ex /2 used for the concentration. Energy behavior vs. 
xpAuj/Ti is similar to the concentration one. To say, a fast 
decrease followed by oscillations. 

Since the in-plane isotropy of the problem, one obtains 
the flow density: i xt = (x/|x|)/ xt , where the scalar func- 
tion I xt takes a similar form to (15): 
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Performing the numerical integrations given by Eqs. (19, 
20) we plot the flow distribution for the above parame- 
ters, as it is shown in Fig. 5. Flow distribution has 
also been normalized by Iq = noPA^/m, being 7 = 
4.57 x 10 18 (cm/s)/cm 2 for the above values of TiAui and 



FIG. 3: Spatio-temporal evolution of concentration, n xt , for 
fT ex /2 = 1. (a) 3D graph near the peak, and (b) contour plot 
showing the line corresponding to the classical velocity. 



jr ex /2. Once again it appears an initial fast decrease, fol- 
lowed by oscillations in dimensionless coordinate. More- 
over, as can be seen, there are oscillations also in time 
and negative values of the flow distribution arise. 



V. CONCLUSIONS 

In summary, we have suggested a new scheme to in- 
vestigate the quantum peculiarities of the single-particle 
dynamics under ultrafast photoionization of a single deep 



FIG. 4: Evolution of energy, £ x t, for the same parameters of 
Fig. 3. 



impurity. Due to the negative contributions to transient 
Wigner distribution of c-band electron, the mean concen- 
tration, energy, and flow demonstrate an oscillatory be- 
havior in contrast to the classical results (see Appendix) . 
We have analyzed the conditions for visible quantum os- 
cillations, when a direct experimental mapping of the 
quantum distribution should be possible. 

Now we turn to the discussion of possibilities for ex- 
perimental verification of the peculiarities obtained. The 
stage of selective single-electron photoexcitation is based 
on the assumption of a low concentration of deep impuri- 
ties: if a bulk concentration less than 10 12 cm -3 remains 
in the near-surface region, one obtains an inter-center 
distance about lfim. The regime of a single-center exci- 
tation can be easily realized with an ultrafast pump fo- 
cused over a submicron scale. Recently, similar measure- 
ments were performed with a single quantum dot [l8[ but 
the photoexcitation into continuum and further evolution 
of distribution was not examined. Perhaps, it is due to 
the complicate problem of the registration of the oscillat- 
ing Wigner distribution. In spite of the sensitive optical 
methods developed recently for optical control of a single 
quantum dot (see [HI, [2(| and Refs. therein) the spa- 
tial resolution remains a complicate task (we use above 
ft/PAui ~ 10 nm). Note that the period of the oscillations 
increases with time as \fi [see Eq. (18) and Fig. 3b] but 
the distribution value (and the response) decreases due 
to spatial spread. Another possibility is to use the scan- 
ning tunneling microscopy [2l| which has nanometer res- 
olution but has to be adapted to time-resolved measure- 
ments with subnanosecond resolution. Note, that we do 



FIG. 5: Evolution of the flow, I xt , for the same parameters 
of Fig. 3. 



not calculate any concrete optical or tunneling response 
supposing that the observed peculiarities will be of the 
same order as the mean values considered in Sec. IV. 

Next, we discuss the assumptions used in our calcu- 
lations. The main approximation is the local time ap- 
proach, so that the edge photoexcitation is beyond of our 
consideration. A more complicate numerical simulation 
is required for this case as well as to take into account the 
Coulomb correlations (excitonic effect). The interaction 
of the electron with the localized hole is essential for a 
near-center region but it should decrease with x. Thus, 
the short-range model used in our calculations of Eqs. 
(9, 11) is enough in order to estimate the photoionization 
decrement. Finally, only the averaged Wigner distribu- 
tion has been considered and a full counting statistics of 
photoionization and subsequent transient propagation of 
electrons [2l| requires a special investigation. 

In closing, a similar theoretical analysis may be de- 
veloped for other cases like photoexcitation of a single 
quantum dot or near-field photoexcitation 0, [l8j] , where 
a similar quantum behavior should take place. We hope 
that these results will stimulate experimental efforts to- 
wards a mapping of quantum peculiarities in the tran- 
sient Wigner distribution. 



APPENDIX A: CLASSICAL EVOLUTION 

This appendix contains the description of the classical 
regime of transient evolution when the generation rate in 
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FIG. 6: Evolution of the classical distributions n x t (a) and 
I xt given by Eqs. (A3) and (A4) for a =6. 



plane scale of excitation, and u>t is the above-introduced 
form-factor. After the integration over q-plane, the dis- 
tribution / PjXt takes the form: 



/ P ,xt oc <5a £ (£ 



-Sex) / 
J — i 



dt'wf exp 



(A2) 

so /p,xt > because of the positive functions un- 
der integral. At t > r ex one obtains the distribu- 
tion as a moving Gaussian peak: / p , x t oc 5j\ s (s — 
e ex ) exp[— (l/2)(x t /Z ex ) 2 ]. The explicit expressions for 

fpix introduced in analogy to Eq. (12) can be written 
through the probability integrals and they have a single- 
peak behavior. 

Restricting ourself to a narrow energy distribution, 
Ae <C Sex, arid taking the integrals over p-plane accord- 
ing the definition (14), one obtains the concentration dis- 
tribution as follows: 



n xt oc exp 



2 G ex ) / 



rt/r ex 



dre 



-2T-* 



l ex \ 1 ex 



(A3) 



where Io(z) is the zero-order Bcssel function of imaginary 
argument. Here we have introduced the dimcnsionlcss 
parameter a = v ex T ex /l ex with v ex — y // 2e ex /m. Within 
the above approximation, the energy distribution is given 
by £ xt ~ e ex n xt . Similar to Eqs. (19, 20) one obtains the 
flow density i xt = (x/|x|)7a; t , where the scalar function 
I xt is written as follows: 



Ixt oc v ex exp 
„-«(*/ 



"2 (z ex ) / 



t/T e 



dre 



*ex \ ' ex 



(A4) 



where h(z) is the first-order Bessel function of imaginary 
argument. 

Performing a simple numerical integration of Eqs. (A3) 
and (A4) one obtains the concentration and flow distribu- 
tions versus the dimcnsionlcss coordinate and time, x/l ex 
and t/r ex , as it is shown in Fig. 6 for a = 6. Since there 
are no oscillation of the classical distribution (A2), the 
mean values appear to be spread monotonically. 



Eq. (8) is approximated by the factorized expression: 

G p .qt oc 5 A s(e - Sex)e- (ql ™ )2 w t . (Al) 

Here 5as(£ — £ex) is the peak energy distribution with 
the half-width Ae placed at the energy s eX i ^ex is the in- 
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